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ABSTRACT 

A relativistic pair-plasma which contains a high excitation of electrostatic turbu- 
lence could produce intense radiation at brightness temperature in excess of 10^'' K by 
stimulated scattering. Important relativistic effects include the broadband frequency 
response of the plasma, and Compton-boosting of the scattered radiation. In radio- 
frequency relativistic plasma, the optical depth can be as small as tens of meters. When 
the plasma wave excitation is one-dimensional and particle distributions have rj_ << Ty , 
the frequency-dependent angular distribution of the emission exhibits characteristics of 
pulsar emission. 

Subject headings: instabilities — plasmas — pulsars: general — radiation mechanisms: 
non-thermal 

1. Introduction 

Radio emission at extremely high brightness temperature is possible by stimulated scattering 
in astrophysical plasmas with a high degree of plasma turbulence. Intensities on the order of 10^'^ 
to lO^'' K are suspected in some astrophysical objects, and an emission process of the form of a 
plasma maser may explain the extraordinary intensity at radio frequencies of pulsars (Melrose 1996) 
and active galactic nuclei (Wagner & Witzel 1996). The following calculation exhibits the intensity 
and spectrum of emission in a relativistic plasma of electrons and positrons under the assumption 
of a uniform and high degree of excitation of plasma electrostatic wavemodes. It shows that the 
path length through the plasma can be relatively small for radiative growth; that Compton-like 
scattering increases the mean emission frequency above the plasma frequency; and the frequency 
bandwidth of the emission is fairly narrow despite the broadband excitation in the plasma. 

The calculation is done in a relativistic regime because pair plasmas derive from extremely 
energetic processes such as gamma-ray annihilation. Simulations of pair plasma creation (Arendt 
& Eilek 2001) show that the pair-plasma distribution functions can be described by a thermal 
parameter which is moderately relativistic. In this paper, the kinetic temperature parameter p = 
mc^ / {kBTjc) is assigned a value of 1/10. 
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Conversion into electromagnetic modes by scattering on plasma waves in nonrelativistic space 
plasma is known to produce radiation at the plasma frequency Up = (47rne^/me)^/^, and wave-wave 
coalescence at twice the plasma frequency (Gurnett et al. 1981). Radiation at these well-defined 
frequencies is due to the narrow frequency response in the plasma. It seems obvious that emission in 
a relativistic plasma is unlikely to follow this plasma emission paradigm. For one thing, electrostatic 
modes in a relativistic plasma exist over a broad range of frequencies, LOp-y/p < lo < ujpj (Godfrey 
et al. 1975a,b; Melrose et al. 2000). Furthermore, scattering by relativistic particles can modify 
frequencies by factors of 47^. Still, the stimulated emission appears to be fairly narrow in frequency, 
but this is a characteristic of maser emission due to frequency-dependent growth acting over many 
growth lengths. 

Although this solution is illustrative of relativistic effects in plasma emission, several assump- 
tions are used to simplify calculation. First, the excitation of the turbulence is characterized by 
a single temperature parameter without regard to specific plasma-streaming or shock-excitation 
mechanisms, or wave cascades. For example, in thermal equilibrium the energy density in plasma 
waves relative to kinetic energy density, E"^ / {STrnkBTx) , is inversely proportional to the number of 
particles in a Debye cube, riA^ = n[c/(ajpy^)]^. This parameter can be quite large for astrophysical 
plasmas. If a radio-frequency plasma (fp ~ 3GHz) manages an equipartition between turbulent 
electrostatic energy and thermal kinetic energy at temperatures of Tk ~ 5 x lO^^K, the charac- 
teristic excitation as described by an enhanced temperature Tnl = {n\\,)TK could be more than 
10^^ K. This coherent enhancement of scattering on plasma waves was recognized by Gailaitis et 
al. (1964), and Colgate et al. (1970). 

Another simplification is to do the calculation without an imposed background magnetic field. 
This justifies the assumption of isotropy in the particle distributions, wave spectra, and wave 
dispersion. Including a magnetic field would tend to make the turbulence one-dimensional in the 
direction of the magnetic field, and would modify the dispersion properties of the electromagnetic 
modes. Thus the plasma maser in a magnetized plasma would not be istropic in direction, but would 
have preferred directions for emission. The effect of making the maser plasma one-dimensional is 
illustrated by an example. 

Finally, the turbulent excitation is assumed to be steady. This energy reservoir in turbulence 
must be maintained by an injection process against the radiation loss, but to model this requires 
an additional kinetic model for the turbulence. The assumption of constant wave excitation can 

be expected to fail when the brightness temperature greatly exceeds the wave temperature, and 
radiation becomes a significant energy sink. Radiative losses will also determine the lifetime of the 
system. 

The calculation proceeds from classical scattering rates between plasma and electromagnetic 
wave modes. The kinetic equation is put into the form of a radiative transfer equation in the next 
section, assuming a kinetic temperature for the plasma and a nonlinear effective temperature for 
the plasma turbulence. The necessary integrals for the scattering coefficients are done numerically 
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as described in section 3. The intense emission which can derive from the maser process is discussed 
in section 4. 



2. Transfer Equation 

We start with a transfer equation between the plasma waves and electromagnetic radiation (see, 
for example, Melrose (1980), section 5.4). For this equation, the electromagnetic wave spectra is 
described by the photon number density function, A'^^(k), and the plasma wave spectra is described 
by the number density function (k) . (The number density spectrum is also called the occupation 
number.) 

By definition, the number of photons with wavevector k in a phase space volume d^k around 
k is given by -/V(k) [(i^A;/(27r)^] . The integral of -/V(k) over phase space will give the number of 
photons per cm^. 

The plasmon occupation number can be formulated from the electric field energy spectrum of 
the turbulent waves, 

J_ f Eir,t)E{r,t) 1 f cLj E*{k,u;)E{k,u;) 

VTJ 47r VTj (27r)3 27r 47r ^ > 

averaged over time T and volume V. Assuming well-defined wavemodes, E{\s.^u) = E(k.)2Tr5[LO — 
a;^(k)], the integration over u gives the following: 



W 



1 E*(k)E{k) 



(2^ 



(2) 



V Att 

Because the term in second brackets is the number density of wavemodes between k and k + dk, it 
is easy to identify the first term in brackets as the electrostatic energy per wavemode. We apply a 
thermal-like excitation of normal modes in the turbulence to derive: 

V 4^ -^ksT^L . (3) 

For electrostatic waves, another degree of freedom is invested in the particle motion. Finally, the 
occupation number is acquired from a semi-classical formula 



1 0, otherwise. 

The inequality describes the frequency range of plasma normal modes in a relativistic plasma. 

With this notation, the transfer equation, in the classical limit, for photons in plasma turbu- 
lence is given by (Melrose (1980), eq. 5.89): 

diV^(k) f d^k' ,3 , , 

/(p) {iv^(kO - iv^(k)} + iv^(k')iv^(k) - k') . 1^1 . (5) 
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The three-dimensional thermal distribution of lepton momentum in a relativistic plasma is given by 
the normalized function f{p)d?p = exp{—p^{p))/Z d^p; here, Z = {4ttK2{p))/ p, K2 is the modified 
Bessel function, 7 = (m^c^)]^/^, and p is the temperature parameter. The spatially-averaged 

number density of particles in the zero-momentum reference frame is n. To write this equation in 
terms of the specific intensity /(w) for a single polarization, use 

I{u)du;dn = nwN^{u/c)c^=^ . (6) 
The resulting equation has three terms. The first term represents spontaneous scattering: 

f = ^^(")(^'«^^^ ' 
where the emission coefficient has been put into the form 



Ai(a;) 



The second term represents absorption 



ds 



where the absorption coefficient is 

/d^k' 
^n/(p)dW(p,k,k') 



(8) 



(9) 
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The final term contributes to stimulated scattering. For a thermal distribution, the momentum 
derivative gives 

df _ c/(p) 

dp 



(11) 



The vector arithmetic can be worked out using the kinematic relationship between frequencies, 

1 -/3 cos 6* 



UJ = LO- 



l-(3cose' ' 

where the angles are between the wavevectors and the electron velocity vector. Thus, 



;i(k - k') 



9/ 

dp 



k^ 



1 - 



UJ 



fiP) 



The stimulated scattering term is 

dl _ Tnl 
""ds T 

The complete transfer equation is 



[A2(u;)-Ai(u;)]/(u;) 



dliiu) iO^ksTNL , 
= Ai , 1- 



ds 



(27r)3c3 



-A2 + %(A2-Ai) 



(12) 



(13) 



(14) 



(15) 
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3. Formulation of the Scattering Coefficients 

The probability for scattering of longitudinal waves into transverse waves by relativistic elec- 
trons is given by (Melrose (1980), eq 4.150; also Melrose (1971)) 

^(27r)3 m^u:'u (1 - • /3)2(l - fe' . /?)2 V ^ \ 

X "(1 - ^ • (l - (y ■ &?) - (1 - /?')(^ -y-^y- ^ , (16) 

where use is made of uj' = k'c. n a7-^^(p, k, k') d^k' has units s~^. In Equation (16), a sum is made 
over polarization in the scattered waves. 

The integration of the scattering rates over d^k' and d^p are done in spherical coordinates in 
which angle 0' measures k' relative to /? and 9 measures /3 relative to k. In terms of angle cosines 
/i = cos 0: 

d^k' = —^duj'dch'du,' ; 

d^p = p'^dpd(j)dn . (17) 

The integral over duj' can be done easily with the delta-function. After integration over d(f)' and 
dcf), the remaining integrals over angle in the emission rate are 

where A and B are algebraic functions of /j,: 

A = 1 + - 4/3/i + ^2 ^ /? V , 

B = 1 - 5/3^ + 2/3^ + - 3^2 ^ 3/3^^ - 2/3V^ • (19) 

The upper cutoff to the dfi' integral is a byproduct of the delta-function, and the high-frequency 
bound on a;' < (1 -|- /3)j(jjp. Thus, 

r 1 _ 

fJ-max = min <^ P 2/37a;p ' (20) 

In order for n'^ax larger than —1, the angles n must be larger than 

r 1 _ 1+13 2-(u)p 

l^min = max I /3 _^ ' (21) 



Similarly, the absorption rate 

A2(a;) 



re^-TcjAnpfip)dpJ^^^^^^-^d,J_^ (13^^^' • (22) 

The integrals in the above equation are completed numerically. Solutions for Ai(a;) and A2(w) 
are shown in Figure 1. The two scattering rates are equal (Ai ~ A2) near the characteristic 
frequency ujp(l/ p^/'^). 
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4. Amplification by Stimulated Emission 

With fixed values for the wave temperature, Tjvl, and the kinetic temperature, Tk, the transfer 
equation is a first order differential equation with constant coefficients which can be solved directly 
as a function of path length: 



I = Ij^e'^^^ir.lc ^ /jvi-^ LsKtiralc _ ^ ^ (23) 



A. 

where the effective stimulated scattering rate is 



Astim(w) = (A2 - Ai)Tnl/Tk - A2 . (24) 

At low frequency (lu ^ Upp^^"^), Ai S> A2 and the intensity is limited to the kinetic tempera- 
ture. Near the characteristic frequency, the intensity saturates at the nonlinear wave temperature. 
At higher frequencies, the intensity increases exponentially. However, the exponentiation rate di- 
minishes at frequencies much higher that the characteristic frequency because Ai and A2 become 
small. 

The high brightness temperatures derived from this theory must be qualified by the fact that 
cooling of the plasma and the turbulence are not taken into account. A full picture of the energy 
balance requires additional kinetic equations for the plasma waves and particle energy distribution. 
However, it is easy to estimate how long the plasma can maintain its relativistic energy against the 
Compton loss. At W^^K, a 10^ cm scale system (such as a pulsar magnetospheric source region) 
will last 10~^s, and a 10^^ cm scale system (encompassing an AGN accretion region) will last for 
lO^s against radiative loss. These simple numbers suggest how maser lifetime might relate to 
microstructure in pulsar radio emission and intraday variability in QSO's. 

To show the magnitude of stimulated emission, the scattering coefficients can be scaled as 
follows: 

^ Zi^ = 1 3 X 10-1° ^ ( ^) ( " ( ^ ) (25) 

c Tk ' Tk \naTC J \ 10^^ cm-^ J \ 1000 cm J ' ^ ' 

Substantial growth in intensity occurs in path lengths smaller than a kilometer assuming turbulence 
temperatures Tnl/Tk > 10^^, as demonstrated in Figure 2. Such a large temperature for the 
turbulence is not implausible from an energy standpoint. The turbulent temperature estimated 
from equipartition between electrostatic and plasma kinetic energy is 

^-^^^^ 6x10-3 nO^^c^'YC^r ■ (20) 



/ c^ \ 



Tk \ujp^/pj V / VP/ 

For a plasma kinetic temperature on the order of 10^'^ K, radiative brightness temperatures greater 
than lO^^K are consistent with stimulated emission from this mechanism. Higher brightness tem- 
peratures are possible with greater turbulent temperature or larger path lengths, depending on 
what limits are imposed by Compton cooling. 
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In summary, the relativistic-plasma Compton maser uses free energy in the form of wave 
turbulence to produce electromagnetic radiation via induced scattering in the rclativistic thermal 
plasma. The high brightness temperature derives from the large value of the plasma parameter, 
n[c/ (wpy^)]^ - generally true in astrophysical plasmas, although less so in laboratory plasmas. The 
maser turbulence conversion mechanism is an alternative to other plasma turbulence conversion 
processes invoking coherent spatial effects or nonlinear waves dynamics (for example, Weatherall 
(1997, 1998); Asseo et al. (1990)), which may not develop due to the turbulence being strongly 
driven or highly inhomogeneous. 

The turbulent conversion process described here might also apply to pulsars. Maser models are 
not new to pulsar radiation physics (for example, Lyutikov et al. (1999); Luo & Melrose (1995)), 
but these masers resemble free-electron masers in which particle beams generate the emission. A 
model invoking the Compton plasma-maser requires further inclusion of anisotropics due to the 
magnetic field, wave dispersion properties in magnetized plasma, and specific mechanisms for wave 
excitation by coUimated particle beams. However, we can simulate these effects by limiting the 
plasma momentum distribution and turbulent wavevectors to a single coordinate axis. 

In one-dimension, the differential scattering cross-sections in terms of the angle cosine /x be- 
tween the emission wavevector k and the z-coordinate axis are given by 

uj^\p,k',k) (27) 

The choice of sign is (+) for electrostatic waves with k' in the negative z-direction, and (— ) for 
waves in the positive z-direction. The calculation assumes waves have relativistic dispersion relation 
uj' = k'c, and the uniform excitation applies to frequencies ujp < uj' < LOpj y^. The one- 
dimensional electron distribution function is given by /(p) dp = exp [—p^ip)] dp/ {2Ki{p)), where 
Ki is the modified Bessel function. 

The transfer equation is the same as before, substituting for the scattering rates Ai 
A^}{dAi/dfl) and A2 — An{dA2/dfl). Here, AQ is the range of solid angles about the turbu- 
lence axis which have the nonlinear excitation. The important difference in one-dimension is the 
directivity of the maser. Figure 3 is a polar plot in angle and frequency. The emissivity in the 
plasma rest frame is largest in directions near perpendicular to the magnetic axis, with some an- 
gular structure due to relativistic thermal velocities. However, the plasma is moving relative to 
the lab/star frame because of the polar cap current flow, and the emission can be expected to be 
beamed relativistically: this is illustrated in the figure with a Lorentz transformation. 

One intriguing consequence of the moving maser is the angular dependence of the peak emission 
for different frequencies. Radio pulsars show a variation in pulse profile (intensity vs. phase) for 
different observing frequency. This is generally interpreted as radius-to-frequency mapping based 
on the presumption that the emission comes from different heights in the polar cap, and that 



dAi{k, p) 

dn 



/• dk' r 
J (2^^ J 



ndpfip) 



l±(3 
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the emission frequency is tied to the local plasma frequency. As clearly shown in Figure 4, the 
maser emission from a single location can produce similar frequency-dependent profiles. Note that 
the spectrum of the emission is found here to have a fairly steep index in the case of uniform 
one-dimensional turbulence, Ii, ~ v^^. 

Finally, we remark that rclativistic temperatures have the effect of supporting wave frequencies 
which are in a broad range about the plasma frequency: still, the highest growth occurs at frequen- 
cies above the plasma frequency. Thus, relativistic effects do not appear to mitigate the puzzle that 
emission at the local plasma frequency in the pulsar polar cap plasma produces frequencies which 
seem high for radio emission (Melrose & Gedalin 2000; Kunzl et al. 1998). 

This work is supported by NSF grants AST-9618408 and AST-9720263. Justin Jayne con- 
tributed to evaluating the numerical integrations. Discussions with Paul Arendt, Jean Eilek, and 
Tim Hankins are gratefully acknowledged. 
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Fig. 1. — Scattering coefficients determined numerically from Eqs. 18 and 22 for isotropic thermal 
particles of temperature UbTk = lOmcP. 
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Fig. 2. — Intensity of plasma maser emission in a electron-positron plasma with kinetic temperature 
of IQmc?. The turbulent excitation is taken to be enhanced by Tj^l/Tk = 10^^. Scaling to a plasma 
frequency of /o = SGHz, the path length corresponds to s = 2 x 10^ cm. 
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Fig. 3. — Contours of constant brightness temperature in frequency space for the Compton-maser 
for one-dimensional turbulence with Tnl/Tk = 10^^, and a path length of s = 1.5 x 10^ cm. The 
frequency is in units of the plasma frequency, /q. The contour levels are 1, 2, 4, and 8 x 10^^ K. 
The contours for both the plasma rest frame, and the lab frame in which the plasma moves with 
7 = 3.8 are shown: the lowest contour is dotted in the former case. Intensity along the three cuts 
of constant frequency are presented in the next figure. 
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Fig. 4. — Brightness temperature vs. angle for three different frequencies for the case of a moving 
plasma. The temperature scale is the same for each plot, with a scale maximum at 1.5 x 10^'' K. 



